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ABSTRACT 

We report on two ASCA observations of the high-mass X-ray binary pulsar 
OAO 1657—415. A short observation near mid-eclipse caught the source in a low- 
intensity state, with a weak continuum and iron emission dominated by the 6.4-kcV 
fluorescent line. A later, longer observation found the source in a high-intensity state 
and covered the uneclipsed through mid-eclipse phases. In the high-intensity state, 
the non-eclipse spectrum has an absorbed continuum component due to scattering by 
material near the pulsar and 80 per cent of the fluorescent iron emission comes from 
less than 19 lt-sec away from the pulsar. We find a dust-scattered X-ray halo whose 
intensity decays through the eclipse. We use this halo to estimate the distance to the 
source as 7.1 ± 1.3 kpc. 
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1 INTRODUCTION 

OAO 1657—415 is one of the most poorly studied eclips- 
ing high mass X-ray binary systems, and yet is potentially 
one of the most interesting. The reason this system is not 
so well known is because it took 15 years for its nature as 
a high mass X-ray binary to be revealed. O AO 1657—415 
was d iscovered with the Copernicus satellite jPolidan et alJ 
I1978T) . A subsequent observation with HEA O-1 revealed a 
pulse period of 3 8 s dWhite fc Pravddll97Slb The 3-30 keV 
Ginga spectrum llKamata et al.ll990f) is a typical X-ray pul- 
sar spectrum with a power law (a ~ 0.6), a high energy 
cutoff (E c ~ 5 keV and Ef ~ 17 keV) and an iron line at 
6.6 keV (equivalent width ~ 240 eV). 

For many years it was unclear if OAO 1657—415 was a 
low-mass or a high-mass system. The mystery was cleared 
up when the CGRO BATSE discovered a ~ 10.4-d binary 
orbit based on long ter m monitoring of the pu lse period 
in the 20-40 keV band IChakrabartv et al]ll993h . A ~ 1.7 
day eclipse of the X-ray source by its companion was also 
seen. This conclusively showed OAO 1657—415 to be a clas- 
sic eclipsing massive X-ray binary system, probably similar 
to Vela X-l, making it the seventh eclipsing X-ray pulsar 
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system to be found. IChakrabartv et all Jl993l) used the or- 
bital parameters to infer that the companion is a supergiant 
of spectral class B0-B6. From the observed value of the spin 
period and its derivative during the spin-up interval they 
deduced an X-ray luminosity of > 1.6 x 10 37 ergs s -1 , and 
thus a distance of > 11 kpc. A more precise X-ray l o cation 
obtained with Chandra allowed IChakrabartv et all (l2002f) 
to make an infrared identification of the companion. The 
infrared properties of the companion are consistent with a 
highly-reddened B supergiant at a distance of 6.4 ± 1.5 kpc, 
implying an X-ray luminosity of 3 x 10 36 ergs" 1 . 

OAO 16 57— 415's position on the Corbet diagram 
JCorbetHl98rJ) suggests that it may not be a typical system. 
The pulse periods of high-mass X-ray binaries (HMXRB) are 
distributed between 0.069 and 835 s, with no evidence for a 
clustering at any particular period. For the Be star systems 
there is a strong correlation between orbital period, P , and 
spin period, P 3 - The supergiant systems have no strong de- 
pendence of orbital period on spin period. OAO 1657— 415's 
38 sec pulsar with its 10.4-d orbital period lies between the 
Be and supergiant systems on the Corbet diagram. The pe- 
riod measurements from RXTE, BATSE, and previ ous ob- 
serva tions give a steady spin-up ti mescale of 125 yr dBavkall 
Il997fi with short-term fluctuations <BavkaJl2()6ch . This sys- 
tem may be in a short-lived phase where it is changing from 
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a wind-accreting system like Vela X-l to a disk-accreting 
system like Cen X-3. Thus, OAO 1657—415 provides an op- 
portunity to test wind and/or disk accretion theory on a 
system which is in a transition between the two states. 



2 THE OBSERVATIONS 

OAO 1657—415 was observed near the middle of its eclipse 
on 1994 March 24, yielding 22 ks of data. This obser- 
vation covered orbital phases —0.001-0.074. It was ob- 
served later (near 1997 September 17) through its eclipse 
ingress for 200 ks. This observation, covering orbital phase 
-0.210-0.011, yielded 67 ks of SIS data. The orbital 
phases in this paper are calculated using the ephemeris of 
Bilds ten et al.l lll997l) . For both of the observations the SIS 
was operated in one-CCD bright mode and the CIS in PH 
mode. The eclipse spectra were significantly different be- 
tween the two observations. We discuss these differences in 
Section mi 



3 RESULTS 

3.1 Timing analysis 

3.1.1 The Light Curve 

Energy-resolved SIS lightcurves and hardness ratios for the 
second observation are shown in Fig. The eclipse transi- 
tion is extended, with the high-energy band becoming fully 
obscured around orbital phase —0.09. The flux in the low 
energy bands decays gradually in the eclipse, suggesting the 
presence of a persistent soft excess. 



3.1.2 The pulse profile 

We derived a barycentric pulse period of 37.39 ±0.01 s from 
the non-eclipse part of the 1997 observation. The pulse pro- 
file from the 1997 observation is qual itatively similar to the 
1-20 keV pulse profile observed by iKamata et al.l (ll990T) 
with the GINGA LAC. Most of the power is in the first 
harmonic. The pulsed fraction also increases with energy, 
as would be expected for a source with such a high absorb- 
ing column density (~ 5 x 10 22 cm -2 ). The dip which oc- 
curs after the main peak (labeled 'A' in Fig.|2l seems to be 
sharper in the present data than in the GINGA data. When 
the 1-10 keV GIS2 and GIS3 data from the first 30 ks of 
the observation are folded into 256 pulse-phase bins, dip A 
has a Gaussian width (a) in pulse phase of 0.0256 ± 0.0054 
at 90 per cent confidence. This corresponds to a duration 
of 0.96 ± 0.20 s. The main dip corresponding to the pulse 
minimum (labeled 'B' in Fig.0 has a width of 0.097 ±0.019 
in pulse phase, or a duration of 3.6 ±0.7 s. Sharp dips in the 
pulse profile have also been seen in Vela X -l fe.g. lChoi et alJ 
I l996f). LMC X-4 (e .g. lLevine et al]ll99lD . and GX 1+4 (e.g. 
iDotani et aflll989h . The pulsed fraction in the pulse profile 
of OAO 1657—415 increases with energy, suggesting that the 
pulsating component of the continuum is subjected to more 
absorption than the non-pulsating component. 
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Figure 1. Top: Energy-resolved SIS lightcurves observed during 
the second ASCA observation. The slow decay of the intensity in 
the low-energy bands through the eclipse is due to the decaying 
dust halo. Bottom: SIS hardness ratios. Note that the lightcurves 
are plotted against orbital phase, while the hardness ratios are 
plotted against time from the start of the observation. Both are 
plotted over the same range to allow direct comparison. 



3.2 Spectroscopy 

The standard screening criteria were applied to the data 
and GIS spectra were extracted from circular regions with 
radii of 6 arcmin. Because this source is near the Galac- 
tic plane, where c ontamination from th e Galactic ridge can 
be significant (e.g. lKovama et alJll989T) . the standard back- 
ground spectra accumulated from blank-sky observations 
could not be used. Thus, we extracted background spec- 
tra from source-free regions of the detectors. To account 
for the azimuthal dependence of the XRT spatial response, 
background spectra for the GIS data were extracted from 
circular regions whose centres were obtained by reflecting 
the source coordinates through the optical axis. To avoid 
oversubtraction of the source, only eclipse data were used 
for the background, and the part of the background region 
less than 8 arcmin from the source was excluded. SIS Bright 
mode spectra were extracted from 4-arcmin circles around 
the source. Background spectra were extracted from the part 
of the chip more than 6 arcmin away from the source, using 
only eclipse data. 
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Figure 2. Top: GIS 2+3 energy-averaged pulse profile. The dips 
in the profile labeled 'A' and 'B' are discussed in the text. Bottom: 
GIS 2+3 energy-resolved pulse profiles. 



3.2.1 The non-eclipse spectrum 

Referring to Fig. the source is essentially uneclipsed 
for the first 30 ks of the observation, or orbital phases 

cf> = —0.210 0.177. We extracted spectra for this interval, 

added the two SIS and GIS detectors, and fitted a model of 
the form 

1(E) = e- r7EN ^(I d E~ a + e - aENH 'I a E' a + I h E- (a+2) + 

simultaneously to the SIS and GIS data. The continuum is 
represented by a power law plus a heavily absorbed power 
law with the same photon index, a. We interpret these as 
representing direct and scattered emission from the pulsar. 
As Thompson scattering is a colourless process the power 
laws have the same index. To model the dust-scattered halo 
we included a power law with a photon index a + 2 and 
the same absorption column density as the direct emission 
power law. We included two Gaussian emission lines, due to 
fluorescence of cold iron and recombination of H- and He- like 
iron, labelled with the subscripts fl and rc, respectively. 

The best-fitting model to the non-eclipse SIS and GIS 
data is shown in Fig. [3] The fit is not formally acceptable, 
with P v (xl) = 6-9 x 10" 3 . Most of the contributions to \ 2 
come from the GIS data. This model fits the SIS data well, 
with P v (xl) = 0-34. 



SIS + GIS Non-eolipse 




S 4 6 8 10 

channel energy (keV) 



Figure 3. Non-eclipse SIS and GIS spectra fitted to the model 
of Equation 



3.2.2 Orbital dependence of the spectrum from the 1997 
Observation 

In order to follow the evolution of the spectrum through 
the eclipse, we divided the data into time intervals based on 
inspection of Fig.Qand attempted to find a consistent model 
that could be used to parameterise the orbital dependence. 
The data were divided into the following five intervals: non- 
eclipse O = -0.210- - 0.177), pre-eclipse (<p = -0.177- - 

0.144), ingress (<j> = —0.144 0.100), early-eclipse (<j) = 

-0.100- - 0.044), and mid-eclipse (<j> = -0.044- + 0.011). 

Fig.^Jshows the SIS spectra at different orbital phases. 
We see the iron line equivalent width increase as the eclipse 
ingress progresses and the continuum decays, while the spec- 
trum becomes harder due to increased absorption. However, 
it is clear that there is a soft component which remains, even 
as the hard part of the spectrum continues to decrease. In 
Fig. 0] the 1994 mid-eclipse spectrum is also plotted. The 
soft component is clearly much weaker in this observation. 

Table Q shows the results of fits to a model of the form 
shown in Equation In these fits we froze the interstel- 
lar column density Nh and the photon index of the direct 
beam from the pulsar a at the pre-eclipse values. Otherwise, 
we were not able to extract meaningful confidence intervals 
from the fits to the early- and mid-eclipse phases as pa- 
rameters became degenerate. As the eclipse progresses, the 
direct and Thomson-scattered emission components are re- 
placed in the best-fitting model by the halo component. We 
thus conclude that the eclipse spectrum is dominated by the 
dust-scattered halo. 

In order to investigate the orbital dependence of the 
iron line complex, we again split up the data but this time 
into 19 finer time intervals. The intervals were chosen to 
provide a compromise between counting statistics and time 
resolution. The exposures were typically 4-6 ks in the early 
parts of the observation, rising to 16 ks for the mid-eclipse 
phase. The lowest number of counts in any spectrum was 
425. The spectra were rebinned so that each bin contained at 
least 30 counts. Each spectrum was fitted to a model which 
had an absorbed power law continuum and the fitting range 
was restricted to the 5-8 keV energy band. Because of the 
degeneracy between the line and continuum parameters, the 
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Figure 4. SIS0+1 spectra at different orbital phases. The orbital 
phases of the spectra correspond, from top to bottom, to those 
in Table^ The mid-eclipse spectrum from the 1994 observation 
is also plotted at the bottom. For clarity, the intensity of each 
successive spectrum has been multiplied by a factor of 0.1. 



absorption column density was fixed for the fits. The power 
law photon index was allowed to vary. 

Fig- 03 shows the results of these fits when the iron line 
complex was modelled by a single Gaussian. When the line 
width was allowed to vary, both the line width and the line 
energy contained considerable scatter due to off-diagonal el- 
ements in the curvature matrix linking them. Thus, the line 
width was fixed at 0.25 keV (the observation-average value) 
for the fits. The line centroid energy increases from 6.45 keV 
in the pre-eclipse phase to 6.55 keV at mid-eclipse. The line 
intensity drops by an order of magnitude as the pulsar is 
eclipsed. To determine the time taken for this transition to 
occur, we fitted a simple continuous model to the intensity 
which assumes that the pre-eclipse and eclipse intensities 
are constant and that the intensity varies linearly during 
the ingress. We find that the transition occurs in 20 ± 9 ks, 
corresponding to a projected distance of 12 ± 6 lt-sec. 

We fitted the iron-line region of the orbitally-resolved 
spectra with two narrow lines at 6.4 and 6.97 keV. As shown 
in Fig. |S| the two lines each undergo a partial eclipse. This 
means that the regions emitting these lines are larger than 
the companion star. 

By looking at how long it takes for the individual iron 
lines to become eclipsed we can estimate the size of the 
emission regions for each of the ionic species (neutral and 
hydrogen- like iron). When we fit our simple intensity model 
to the line intensities in Fig.[§]we find that the line intensities 
decrease by factors of 4.5 ± 0.6 and 2.2 ± 1.2, in intervals 
of 30 ± 19 and > 23 ks, for the 6.4 and 6.97 keV lines, 
respectively. These time intervals correspond to projected 
distances of 19 ± 12, and > 14 lt-sec, respectively. Thus, the 
6.4 keV line undergoes a deeper and more abrupt eclipse 
than the 6.97 keV lines. 80 per cent of the 6.4 keV emission 
is located within 19 lt-sec of the neutron star. 
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Figure 5. Orbital dependence of the iron line energy and in- 
tensity. The power law photon index is a, Ip^ is the power law 
intensity in cts s — 1 cm -2 at 1 keV, E is the Gaussian line centroid 
energy, and / is the line intensity in cts cm -2 s — 1 . 
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Figure 6. Orbital dependence of the iron line intensities for two 
narrow lines at 6.4 and 6.97 keV. The power law photon index is a, 
7pi is the continuum power law intensity in cts cm -2 s — 1 keV -1 
at 1 keV, and the line intensities are in cts cm -2 s — 1 . 



3.2.3 The mid-eclipse spectrum 

For this orbital phase we have data from both the 1997 and 
1994 observations. We therefore discuss these together in 
this section. For the 1994 observation the intensity in the 
2-10 keV band is 6.5 x 10~ 12 erg cm -2 s _1 . The iron line's 
equivalent width is 5.68 keV. The continuum is weak but 
a single power-law continuum is not acceptable. A model 
which gives an acceptable fit has a continuum composed of 
an absorbed power-law and a second power-law whose pho- 
ton index is constrained to be equal to that of the first one 
plus 2 with less absorption, and two Gaussian emission lines 
whose centroid energy and width are free (see Fig.|7J. Note 
that this model has a simpler continuum than the model 
described by Equation and may be expressed as 



1(E) 



' N ">(I h E-<- a+2) + e -"EN HlLE - a 

+ Ef =1 ^-< B -^ 2 /^ 2 ) 



(2) 



ASCA Observations of OAO 1657—415 and its Dust- Scattered X-Ray Halo 5 



Table 1. Best-fitting model parameters for fits to the orbital phase-resolved GIS and SIS data. The model 
is of the form shown in Equation Q The interstellar absorption column and the photon index of the direct 
power law are fixed at the pre-eclipse values (Njj = 7.61 ± 0.35 cm -2 and a = 1.361 ± 0.043). Note how the 
direct and scattered continuum components (I d and I a ) are replaced by the dust-scattered halo component 
(Ih) as the eclipse progresses. 



Parameter 




non-eclipse 


pre-eclipse 


ingress 


early-eclipse 


mid-eclipse 


SIS Exposure (ks) 




20.8 


24.8 


16.0 


40.5 


29.0 


SIS count rate (cts 




3.30 ±0.01 


1.62 ±0.01 


0.451 ±0.005 


0.205 ±0.002 


0.142 ±0.002 


I d (10~ 3 cts cm -2 


s" 1 at 1 keV) 


33.5 ± 2.2 


8.8 ± 1.4 


2.90 ±0.61 


0.39 ±0.35 


< 47 


N m era- 2 




14.13 ± 0.71 


30.0 ± 1.8 


61.3 ±6.0 


87 ±54 


35 ± 12 


I s (10~ 3 cts cm~ 2 


s" 1 at 1 keV) 


92.9 ± 3.2 


95.2 ±3.3 


30.9 ±3.1 


2.3 ±1.3 


2.07 ±0.48 


Ih (10 -3 cts cm -2 


s" 1 at 1 keV) 


< 13.1 


53.3 ± 1.0 


52.3 ±5.6 


35.8 ± 3.5 


56.0 ±2.8 


E a (keV) 




6.398 ± 0.021 


6.356 ± 0.056 


6.361 ±0.053 


6.44 ±0.47 


6.366 ± 0.050 


CTfl (eV) 




195 ± 54 


227 ± 88 


< 139 


< 131 


< 149 


7 fl (10- 4 ctss" 1 ) 




15.1 ± 2.3 


10.1 ± 5.3 


2.54 ±0.63 


1.37 ±0.48 


1.94 ±0.90 


EW a (eV) 














E IC (keV) 




7.15 ± 0.14 


7.05 ±0.45 


6.638 ± 0.038 


6.87 ±0.19 


6.76 ±0.16 


(T rc (eV) 




< 322 


< 0.350 


214 ± 135 


233 ± 116 


207 ± 77 


irc(10- 4 CtS S" 1 ) 




3.1 ± 1.8 


< 8.95 


3.80 ±0.78 


1.04 ±0.52 


1.07 ±0.73 


EW VC (eV) 














X 2 Jd.o.f 




1.210/367 


1.106/355 


1.046/264 


1.275/248 


0.884/174 






3.52 x 10~ 3 


0.0826 


0.292 


2.17 x 10~ 3 


0.862 



for the case where there are N Gaussian emission lines. Here 
Nhi is the interstellar column density and is the inten- 
sity of the power law due to the dust halo. Nhi is the local 
column density for the absorbed power law, due to scatter- 
ing of the direct pulsar emission in the stellar wind. The 
absorbed power law's intensity is I s - The photon index of 
the absorbed, harder power law, a, was fixed at 1.4 as it 
was found to be unconstrained by the data. The best-fitting 
parameters are shown in Table [5] The 1994 observation ap- 
pears to have caught the source in a low-intensity state. 

In contrast, the 1997 mid-eclipse spectrum has rela- 
tively strong continuum emission. The 2-10 keV flux is 
1.1 x 10 — 11 erg cm~ 2 s — l . Again, a model with two Gaus- 
sians was required to adequately fit the iron K-line region. 
The results are shown in Fig. □ and Table H 



3.3 Spatial analysis 

From the slow decay of the soft-band lightcurves through the 
eclipse we suspected that an extended dust-scattered halo 
might be present. In order to investigate this we extracted 
SIS and GIS images in various energy bands. Surface bright- 
ness profiles were generated for the SIS and GIS data. The 
SIS was operated in one-CCD mode which is not optimal for 
studying structures extended over more than a few arcmin. 
This means that if we try to integrate a surface brightness 
profile over more than about 4 arcmin, we start including ar- 
eas off the chip. Combined with the considerable azimuthal 
and energy dependence of the XRT, this complicates spatial 
analysis. Thus we used the ASCA-ANL software to gen- 
erate ray-traced simulated images which we could compare 
with the data. In each case we simulated the image expected 
from a point source at the centroid coordinates of the ob- 
served source, and with the same spectrum. We then inte- 
grated the data to produce a radial surface brightness profile 
(SBP). We integrated the the ray-traced image in exactly the 
same way, from exactly the same region, to produce a SBP 
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Figure 7. Top: SIS0±1 mid-eclipse spectrum and model for 1994 
observation with fit residuals. Bottom: SIS0+1 mid-eclipse spec- 
trum and model for 1997 observation with fit residuals. The data 
are for approximately the same orbital phases (—0.001 — h 0.074 
for the 1994 observation and -0.044-+0.011 for the 1997 observa- 
tion) and the same spectral binning has been applied. The model 
parameters are as shown in Table |21 Note the much stronger con- 
tinuum in the 1997 observation. 
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Table 2. Best-fitting model parameters for the mid-eclipse data of the 1994 and 1997 observations. The 
model is of the form shown in Equation|2]with two Gaussians representing iron fluorescence and recombina- 
tion emission lines. The data are for approximately the same orbital phases (—0.001 — h 0.074 for the 1994 
observation and —0.044 — h 0.011 for the 1997 observation) and the same spectral binning has been applied. 
Errors are for 90 per cent confidence in one interesting parameter. 



Parameter 




1994 Value 


1997 Value 


SIS Exposure 


(ks) 


29.0 


40.9 


SIS count rate 


(cts s" 1 ) 


0.142 ± 0.002 


0.026 ± 0.001 


N m 


(10 22 cm" 2 ) 


10.0 ±2.6 


6.67 ±0.66 


Ih 


( cts cm -2 s — 1 at 1 kcV) 


(6.1 ± 2.0) x 10" 3 


0.028 ±0.012 


Nm 


(10 22 cm" 2 ) 


62 ± 16 


< 6.20 


Is 


( cts cm -2 s _1 at 1 kcV) 


(4.4 ± 1.5) x 10~ 3 


(1.57 ±0.23) x 10 


a. 




1.4 (fixed) 


1.4 (fixed) 


Efi 


(kcV) 


6.416 ± 0.018 


6.474 ± 0.046 




(eV) 


41 ± 30 


< 118 


h 


(10~ 4 cts cm~ 2 s" 1 ) 


1.54 ±0.27 


1.96 ±0.96 


Ere 


(keV) 


7.13 ± 0.15 


6.83 ±0.17 


Ore 


(eV) 


260 ± 190 


< 380 




(10 -4 cts cm" 2 s" 1 ) 


0.92 ± 0.39 


2.3 ± 1.1 


xS/d.o.f 




0.983/45 


1.074/46 



which represented the point spread function (PSF) expected 
from a point source at the same position with the same 
spectrum. We parameterised the PSF with three exponen- 
tials and fitted for the widths and relative normalizations. 
We then fitted the resulting three-exponential model, along 
with a constant background, to the SBP obtained from the 
data, keeping the widths and relative normalizations fixed. 
The fits were performed only for the innermost and outer- 
most one arcmin of the SBP, where the point source and the 
background should dominate, respectively. 

Fig. |H| shows surface brightness profiles from SIS0 for 
the 0.5-3 and 5-10 keV energy bands fitted to the PSF and 
a constant background. Comparing these, it is clear that the 
low-energy data require a spatially extended source, while 
the high-energy data are consistent with a point source and 
the XRT PSF. We interpret this as evidence for a dust- 
scattered halo. 

We can obtain a rough measure of the halo intensity 
from the residuals of the fits but for a more quantitative 
analysis we decided to fit the SBPs to model haloes, con- 
volved with the ray-traced PSF. We fitted the SBPs to a cen- 
tral point source, constant background, and dust-scattered 
halo. In the Rayleigh-Gans approximation the radial depen- 
dence of the model halo is [ji(x)/x] 2 , where ji(x) is the 
spherical B esscl function of order 1 and x — (4na/X) sin(t9/2) 
( Havakawal ll970l) . Here a is the grain radius, A is the X-ray 
wavelength, and is the angle through which the dust grain 
scatters the X-ray photon. The Rayleigh-Gans approxima- 
tion is valid provided that 4-7ra/A sin(t9/2) <C 1, which means 
in practice that we should restrict our analysis to energies 
above 1 keV (e.g. lSmith fc Dweckll 998). The observed angu- 
lar distance of the detected X-ray photon from the centre of 
the image is a — (l — z)0, where z is the distance from the ob- 
server to the dust grain divided by the distance from the ob- 
server to the source. The model was convolved with the ray- 
traced PSF for the fit. To simplify computation, instead of 
using the Bessel- function form above we used the Gaussian 
approximation of lMauche fc Gorensteinl l)l986l) for the halo 
shape. Figure [5] shows typical fit results for an SIS0 eclipse 
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Figure 8. SIS0 surface brightness profile for the mid-eclipse phase 
of the 1997 observation. The data are fitted to the XRT PSF 
plus a constant for radii less than 1 arcmin and between 6 and 
7 arcmin. Top: 0.5-3 keV. Bottom: 5-10 keV. 
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Figure 9. 3.0-4.5 keV surface brightness profile for OAO 1657- 
415 for the mid-eclipse phase of the 1997 observation. The data 
have been fitted to a central point source, constant background, 
and dust-scattered halo. 



halo. We divided the SISO and SIS1 eclipse data into energy 
bins ( 1.5-2.5, 2.5-3.5, 3.5-4.5, 4.5-5.5, and 5.5-7.0 keV) and 
obtained the best-fitting value of 4-7ra/A for each energy bin, 
assuming that the dust grains were distributed uniformly 
along the line of sight. We then plotted 4-7ra/A against en- 
ergy and the slope yielded the value 0.095 ± 0.005 fim for 
the size of the dust grains. 

The size distribution of the grains may be found either 
from the time d ependence of the integrated halo spectrum 
\Xu et alJll986T) or from the surface brigh tness profile at a 
given energy iMauche fc Gorensteinlll986l) . Since we find it 
easier to separate the central source from the halo in dis- 
crete energy bands by deconvolving the surface brightness 
profile from the point spread function, we adopt the latter 
approach. We do not have good enough counting statistics 
for the former method. 

We split the SISO and SIS1 eclipse data into 20 ks 
(gross) bins and extracted surface brightness profiles for dif- 
ferent energy ranges ( 1.5-2.5, 2.5-3.5, 3.5-4.5, 4.5-5.5, and 
5.5-7.0 keV). We obtained the halo intensity for each time 
and energy bin. We then fitted these halo intensities with an 
exponential function for each energy bin. The spatial reso- 
lution of the GIS is such that the profiles are more smeared 
out. This means that the GIS values for the decay time con- 
stant have larger error bars and in some cases are not con- 
strained by the data. Thus, although the GIS results for the 
decay time constant are consistent with the SIS results, we 
use only the SIS data to obtain the decay time constant. If 
we assume the dust is uniformly distributed along the line of 
sight, we expect the halo to decay exponentially with a time 
constant of td{E) = 13.6-E -2 D8ao i d, where E is the photon 
energy in keV, D$ is the distance to the source in units of 
8 kpc, and ao.i is the dust grain radius in units of 0.1 /im. A 
fit of this form to the data is shown in Fig. 1101 This results in 
a distance of (7.9±0.9)oo.i kpc for OAO 1657-415. We have 
implicitly assumed here that the eclipse is a step-function. 
The eclipse ingress is much shorter than the eclipse dura- 
tion and, as our measurements of the source intensity start 
just before eclipse ingress, it seems reasonable to assume a 
constant intensity prior to the eclipse ingress. However, this 
estimate of the distance is sensitive to the distribution of 
the dust along the line of sight. If the dust is all in a sheet 
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Figure 10. Halo decay time-constants plotted against the inverse 
square of the energy. The SISO, SIS1, GIS2, and GIS3 data are 
represented by open circles, filled circles, open triangles, and filled 
triangles, respectively. The data have been fitted to a function of 
the form t^E) = tE~ 2 and the best-fitting value is r = 18.6 ± 
2.3 d. 



half-way to the source our estimate for the distance becomes 
(12.7 ± 1.5)og. x kpc, while it is (4.56 ± 0.54)a§.i kpc if the 
sheet is at either 10 or 90 per cent of the distance to the 



4 DISCUSSION 

4.1 Intensity States 

We have observed OAO 1657—415 during the latter half of 
its eclipse in two distinct intensity states and have found dif- 
ferences in the spectra. The low-state spectrum is consistent 
with an absorbed power law with photon index 0.36 ± 0.19 
due to scattered and absorbed emission from the pulsar, 
combined with a dust-scattered halo and fluorescent emis- 
sion from cold iron. The high-state spectrum has a stronger 
continuum which is due to scattering of the direct pul- 
sar emission in the stellar wind combined with the dust- 
scattered halo. In the high state the recombination line due 
to ionised iron is stronger than the fluorescent line. In the 
low state the fluorescent line is stronger. This difference be- 
tween the two states is most likely due to increased pho- 
toionization of the stellar wind in the high-intensity state. 

4.2 The spectra 

The aim of the 1997 observation was to use the eclipse of 
OAO 1657—415 to study the outer atmosphere of the OB 
companion star and its vicinity by mapping out the X-ray 
emission in the system. As the eclipse ingress progresses, 
the X-ray emission must propagate through the wind and 
outer atmosphere of the OB star. From the progress of the 
partial eclipses of the iron emission lines we conclude that 
the 6.4-keV fluorescent emission originates near the neutron 
star while the 6.9-keV recombination emission comes from 
the extended wind of the companion. 

We find that a single power law cannot describe the con- 
tinuum spectrum outside of eclipse. We obtained an accept- 
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able fit using two power laws with the same photon index, 
one of which has a larger absorption column density than the 
other. This model is equivalent to a partially-covered power 
law. Similar two-component cont i nua have been observed 



in Cen X-3 jEbisawa et alJll996t |Audlevlll997j), Vela X- 
1 <Sako et alJll999l) . and~GXl+4 jEnddl200lD . ISako et alJ 



(1999) claim that the more absorbed component is the direct 
beam which is absorbed by the stellar wind surrounding the 
neutron star and that the less- absorbed component is due 
to scattering in the stellar wind. lEbisawa et al] il996T> claim 
the opposite: that the direct beam is the less-absorbed com- 
ponent. One way to distinguish between these two pictures is 
by pulse-phase-resolved spectroscopy of the continuum. We 
would expect the component that is due to the direct beam 
to pulsate at the pulsar frequency. As long as the light travel 
time over the extent of the stellar wind is at least a signifi- 
cant fraction of the pulse period, p ulsations in th e scattered 
component should be smeared out. Endo (2001) did this in 
the case of GX1+4 and fo und that both components origi- 
nate near the neutron star. lAudlevI (119971) came to the same 
conclusion for Cen X-3. When we fit the pulse-phase resolved 
non-eclipse SIS data with the two-power law model, we find 
that the intensities of the power laws have essentially the 
same dependence. We therefore conclude that both of these 
components must originate at the same place, significantly 
less than 38 lt-sec away from the neutron star. 

4.3 The dust-scattered halo 

Scattering of X-rays by interstellar dust leads to the forma- 
tion of a halo whose intensity, shape, and spectral properties 
depend on the distance to the source, its intrinsic spectrum, 
the d istribution of du s t grains, and their size and com posi- 
tion jHavakawalll970D . lTriimper fc Schonfelderl <1975ft pro- 
posed that the time-dependence of a dust-scattered halo 
from a variable X-ray source could be used to determine 
the distance to t he source and th e distribution of dust along 
the line of sight. IXu et al] il986T> suggested that high-mass 
X-ray binaries with their abrupt total eclipses that occupy a 
large fraction of the orbit would be good candidates for this 
kind of study. iMitsuda et al] dl990l) used lunar occultation 
observations with Ginga of two X-ray sources to study the 
energy dependence of the dust halo and inferred the pres - 
ence of iron in the sca ttering grains . |Pav fc Tennantl (1991) 
applied the model of iMolnar fc Mauchd lll986t) for the X- 
ray halo of a time-varying source to the low-energy eclipse 
lightcurve of Cen X-3 and obtain ed a distan c e of 5 kpc 
and a dust grain radius of 0.33 /im. IWoo et al] (119941) used 
the imaging capability of ASCA's SIS to measure the sur- 
face brightness profile of the dust-scattered halo in Cen X- 
3 and infer some properties of the dust grains. Applying 
the Rayleigh-Gans approximation, from the lower than ex- 
pected relative intensity of the halo they concluded that 
the dust grains were loose aggregates ra ther than solid par- 
ticles. However, ISmith fc Dweckl (I1998T) later pointed out 
that the Rayleigh-Gans approximation for the differential 
scattering cross-section overestimates the halo intensity be- 
low 1 keV and that when the more exact Mie solution is 
used the data may be consistent with solid grains. Later, 
iDraind i2003l) pointed out that the median scattering an- 
gle for the dust halo would be small enough that most of 
the halo photons would be seen in the core of the ASCA 



image which is dominated by the central point source. This 
would have led I Woo et al] ]l994l) to underestimate the halo 
intensity and is likely to be a problem with t he cu rrent ob- 
servation of OAO 1657-415.. IPredehl et all ]2000h demon- 
strated that Chandra could be used to determine the dis- 
tance to Cyg X-3 to within 20 per cent by correlating the 
lightcurves at different radii from the source. Cyg X-3 only 
undergoes a partial eclipse while high mass X-ray binaries 
such as OAO 1657—415 undergo total eclipses for large frac- 
tions of their orbital periods. This makes a high-mass X-ray 
binary such as OAO 1657—415 a promising candidate for 
this kind of study. 

We have found a decaying dust-scattered halo in the 
eclipse of OAO 1657—415 (see Figure^. The half-power ra- 
dius of the halo is about 4 arcmin, which is not much larger 
than the XRT point spread function that smears out the halo 
and mixes in flux from the core (see Figure [§J. The X-ray 
spectrum of OAO 1657—415 in eclipse has three components: 
the dust halo, X-rays from the eclipsed pulsar scattered by 
the primary's atmosphere, and emission (mostly iron-line re- 
combination) from the primary's atmosphere. Only the dust 
halo will appear to be extended in an X-ray image. Because 
this source is heavily absorbed (Nh ~ 5 x 10 22 cm -2 ) the 
hal o is mostly seen abo v e the ROSAT pass-band. 

IChakrabartv et all dl993fl used the observed values of 
the pulse period and its derivative during a steady spin-up 
interval to place a l ower limit of ~ 11 kpc on the distance 
to OAO 1657-415. IChakrabartv et al] (120021) found a dis- 
tance of 6.4 ± 1.5 kpc from the reddening of the optical 
counterpart. The distance we estimated from the decay of 
the halo, (7.9 ± 0.9)ao.i kpc, lies between these two values. 
. It is interesting to compare the case of Cen X-3 where 
iDav fc Tennant) dl99l]) obtained a distance of 5.4 ± 0.3 kpc 
which is a little lower than the lower limit derived from the 
relation betw een temperature and luminosity for the BO 
star (6.2 kpc: lKrzemiriskilll974l) . The value we obtained for 
the maximum dust-grain radius is 0. 095 ±0.005 )im which is 
smaller than the value obtained by IDav fc Tennant] (11991]) 
(0.34 ± 0.01 /im). There may be systematic effects resulting 
from our simplistic assumptions which causes us to under- 
estimate the distance and overestimate the grain size. We 
would like to point out, however, that the distance and grain 
size estimates are independent; the measured halo intensity 
does not depend on the assumed grain distribution. Our es- 
timate of the distance is sensitive to the exact distribution 
of the dust along the line of sight; we have assumed that 
it is uniformly distributed. Our results also depend on the 
assumption that the intrinsic intensity and spectrum of the 
source were constant before eclipse. 



5 CONCLUSIONS 

The eclipse spectrum of the source in its high-intensity state 
is dominated by the decaying dust-scattered halo. We esti- 
mate the radius of the grains to be 0.095 ± 0.005 fim and 
the distance to the source to be (7.9 ± 0.9)a§ x kpc, where 
ao.i is the grain radius in units of 0.1 ^m. Using our esti- 
mate for the grain size, the distance to the source becomes 
7.1 ± 1.3 kpc. We find that the non-eclipse high-state con- 
tinuum spectrum cannot be modeled by a single power law. 
A more-absorbed component due to scattering by material 
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near the pulsar is required. The iron lines underwent par- 
tial eclipses during the 1997 observation. We estimate that 
80 per cent of the 6.4-keV fluorescent emission originated 
less than 19 ± 12 lt-sec away from the pulsar while the 6.9- 
keV recombination line was mostly emitted by the extended, 
ionised stellar wind. During the 1994 observation, the con- 
tinuum emission was much lower and the iron emission dur- 
ing eclipse was dominated by the 6.4-keV fiourescent line, 
which suggests that the source was in a low intensity state. 
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